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de Intredyetien. It 18 well knewn in the theory of 
Piane atrain in orthetroepic anteriale” that if F satiszies 
the equation 


We r/ex*) + aya (o4ze2%y?) oa(2*F/oy") =0 a) 


where d, Bo , and U depend upen the eiaatio constante of 
| the material, then the atresses are given by 


f.. = B2p Lay? Tyy = 27h x2, Tay = ee o28f xy. (2) 


The solution aust also satisfy the boundary congitions 
at the surface. These will be called the surface boundary 
eend) tions. Fer laminated gaterizl 1% im also required tha% 
the normal and shear stresses, and the displacements, be eon-« 
tinugus at every puint on the aurfaces separating the layers. 
Those will be culled internal boundary atregs conditions ant 
Anterye) boundery displacement conditgons respectively. . 


fhe strain componente are defined by 


@yy 7 2U/2K, yy = IV/2Y, Oyy = Hav/dy + Bv/dz). (3) 


where u and v are the displacements in the x ana y Alreetiang 
respectively. 

Using E, for the Young's moduli, Wy for the Poisson's 
Katies, and yy for the moduli of rigidity, H2o8t lav con 
be written 


Coy = WT, (4) 


whore. | 
oe he Tx ex) By 
6 2 (Cgy * Cggcgy)/Bz = (yz + ya ply 
| Y = (1h - yagny)/By 
V = V/A. (5) 


Soaving..( 4) for Tyy, yy and Ty givea 


where 
f= p/lar -/*) 
td = nl» 2. % 
= Fiery -B*). (7) 


Let Hs i¥/ol#y aad x= (v~a}/tewi®, uben equation (1) 
gan be @ritten 
ae/az* * ant 342/3x%3,7) + SP/a7,* = Os (8) 
This equatien is satisfied by 


P= RIP, (ty) * Polda)) (9) 


where ; 
2, =x iiiy = ze alk+ (K? . 1) yt, 
Bo = X + ibby = x > ilk - (x2 . 1)3] b> {10) 
with 
6, = (v/a) BLK + (x? - ee 
be = (V/a) 4K - (nx? - 2838, (23) 


If wo let R, denste the resi part of a*7,/a2,*, and & the 
imaginary part, we have 


tx = - 85R, ~ 58k, 

Typ. PE eee 

txy = 83, + bg (32) 
ant 

xx = - 1) - “gRo 

®yy = ea 8 + SR 

Syy = vé)8y + Vboko (15) 

2 e: 20¢ 

where €, = aby + B and $i = Aby + Yr. (14) 


§ $4 + Ss, zs 2¥ i zab $ op6e * Y¥ . 252 f 


ne have 4 S Sctav> 6,35. | (15) 


end &1480 
$252 = Y¥/a. (26) 


52 + 62 = av - ae. (17) 


In this paper #6 consider rectangles consisting of three 
layers cemented together. The two suter layers, or faces, are 
similar end of thickneas (b-e). The thistnoss of the innex 
leyer, or core, is 2c. Tne width of each layer is 2a. The 
pianes of s@®paration of the layers are at y = + a. Since the 
problem is one of plane strain, the length of the body in the 
2 direction is considered to be infinite. We put bens and e=-we. 

The internal boundary stresa condititns are thet Toy and 
Fo3 58 somtinuous at y = + ¢. Ths internal boundary Maviece 
ment conditions are that u and v be continuous at y = + ¢c. 
in part 2 we gee that these internal boundary conditions 
6liminate many solutione of egimiien (8). 

In part %3 we sive a method of obtaining approximste 
aclutions which satiafy surface boundary esnditions on the 


Streases, oF surface boundary eanditions tn =} 


a oie = - Sas sie Bmenti 5 & = _ % Mtoe 

BiaGew~ente ehich Sit given ave specified t2 be fers. This 
we —s - ~ ~e . Ra ms. 

method is based on ike ecinmnert sf fumctisn apess. 


In part 4 the saqustiscsa ef silastic; 
finite squationa t5 xive a metnod af obtaining uppresimate 


f@olutione which #atisiy surface Doundary conéitians on the 


Stresses. 


ea 


Yee een! ne 
a= re rere rrenmerrnccmnenre 


For an ex@uple we take the saterial to be apruce with the 
longitudinal azia in the direction of the x-axis in the middle 
iayer and the transverse axis in the direstion of the x-axie in 
the outer leyers, tue radial axis being in the direetion of the 
y-sais in both layers. We take the constante ef the metsrial 


te be aa follswe: 


B= 1.95°20° 1b/in“ Ep = 0.128°10° 1p/in? —-& = 0.0691-10° ip/ie? 
‘eds eaae® va Avy = 0.0720:108 ib/in? An = 9.0045-16 1b/an2 , 
Tra = 9.450 Tar = @.559 Ty = 0.0194 

ey, = 9.0360 Vor = 0.501 Tip = 0.539 (28} 


Using the valuea given in (18) #e 


<, D 9 
~ be 
m \ is 


~ 
ui 


0.607°107° in#/ib 


= 0.314°10°9 in@/1p 


6.48°107° 4n2/1p 


4.81°107© an@/ib 


th 


2.00: 1N9O en2@/op. 


0.696°1075 in®/ap 


11.8°10°° in?/ibd 


6.72°107© in@/ib 


Ay 


Bz 


ns 


eaicgulates 


u 


14.3°2075 in@/ib 


4,48°10°© in2/ib 


= 7.69°107© in@/ip 


= 108.7°107© in®/ib 


212°10°© in2/1d 


§.01°10°© in®/id 


72.8°107© in®/ip 


= 7,85°107° in@/1ipb 
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2. Kyact Polynomial Solutiong. Let 
nee 
od eas othe * IM )2) 
Fo = £iGn * iNa)zg””*. (2) 


n, indicating the layer by meana of a second audecript, with 
1 for the core and 2 for the faces, the internal ooundars 


atreaa conditions can be eritten 


<a-& | 5 
“Eo (-2) BPC ae pe2) (ewer) (“EP)eP(51 4 PP ee 1 * $0)PQgep 1 


Pp even 
~ §29¥Payn 2 : 6208 cep, 2) * (2) 


me as = -1) 8P-8( gepo2) (eps) ( cP (51 Peep, * 62," * ace, 1 
p rae 
~ 6)2"" Pen 2 - Sa0"*4ae.5 2) = 0 (3) 


. “En -1) 89-8, erpe2)(iceper) (PP) eP(6 4. Men a + S22 Mery a 
p 08d 


- 51 27%leep 2 - S20PNuvp,2) = 0 (4) 
=a;* BP KeP) Pc 5, Ply 9, PPD 
PEO (-1) °° (urpe2) (eeprh) tg” oP 544" “Beep 2 * 627" “Nigep 1 
p oven 


4 


wine wad 
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Tne internal boundary dieplasement conmlitions san be written 


Bek 
£ (-2) 8P( ce pee) (P74) oP (635 P ea Peep 1 + 55, Pea,Qnen 3 
Pp even 4 4 


~& 2° 2  aPaep,2 > S22” 22@kep,2) = 9 (6) 


1 gtpr) 
MES (BEPC epee) (POP (5, yP Sy aPee ps * S21 7S2iduep,1 
p even 


an 
ba 
= 


~ S59? SpoPeep 2 - San” Sga@q+p, 2) = 0 


EC -2)#P Bi pe2)0P( 55, PPPeaaepy + BoxP eoiay 1 
- ~ 5y_P*4es oP, a - Saqh*2en08, 2) = 0 = 
“pet at k+pr2)( gi )oP(8)1 P63 Mee at 821 °€o Mae, 1 

~ 5:2” €,olgen 2 - SaoP€goNeep 2) = 0 (9) 
aman =) BP-*, arpr2) ( K*ptl) P| 613 Sy Beep, + 63,"5a,R +p, 1 
p oé 


- 6,275) aMep,2 - Se2PS2aNe+p,2) = 0 (RO) 


+ BP-8( 040) 5P(5._P 2 P 
pe CR br 2)97 (874 °S a Mp1 * 81 Saag,1 > Sia Sill, 2 
aa 


im equations (2)-{7} the range of K ig given by O < E <a, 
in equation (9) it is given by -1 < & < m-1, and in equation 


(10) by 1 < & < mri. 


Equations (1) - (11) have ne non-trivial solutions form > 2, 


For g = 2, we get 


Poy = Q,5 = Ms = Nyy = 9, fori =i, 2; Ji, 2. Si 

2 2 os 

‘os b22 ~ 42 bs “21 - 62| S02 : 
62, - 832 {42 $a. - Yat A te 


¥ ve ee ldo E11 - fal Ano . 

S01 ~~ e2 = 2 2 = ¥ | A a 

991 7 911 [42 $12 2 5 
a Eno - Q Cig a 8 
foo 4 > | 22 - P| 2 ee i 12 7 Ay 
PL $22 = Ny A Py $ 12 °* Vy 

8 € a €12 592 € ~ & $22 * 

ma =" ze see Beg, Mai? F 2 Nog, Mag S «Moa 

cs) aaa © 5 el “as * “id 12 : 


Noo 


log = - £22 tog - 202 $22 S22 SaPA-A _ Barter) 


S12 $42 “4 - “2 


S22 €22 ~ ©)2 
OL ge ea ee, oe 44529) 
S24 


€2 & a 
o22 <22 “ (Noo + KaNoo) 


g 
il 
é 
| 
fd 
ih 
fd 
bt 
é 
fy 
~~ 
ps 


x, = 202(2-Ao2eun _ SpenSan BAK _ Brena , See 
4 Ay Uyp~% ay Ag Y-4 


3s 


10 
Eo2~",2 fs Po-Ale €207€12 aoe 
Aggy 3 ee ago S02 + 24 § EBS Noox ie 
€22-Q,2 4, €1-Fo1 Qo2 €a27€y2 “3 
fxx2 7 2 7 ae 820 OH RaoXd oe 
2 la Suny 4 ' age 
fey = Type = - 2 eget [E92 | Soe LA 
n re 
E95 E | Br, -2y a & a 
_ gaae £22792 (As ~ ee py ae 8 
Tey * 48° S35 “He, ay ting ~ 1289 “ay 2ay® ‘ 
“gas _ a ak. 2a oy 
Tzy2 = 40°S20 (act ap Nee - 12692 nag taal” © 
Por the given numerisal example: 2 
= 2 ie a 2 S 
Bo» =~ = 1.18Noo - 120e Noo: Roy = 5-59Ngo * 22. 4e8 Woo 
Myo = = @-0505Ngp - 5: o8e°N.,., My, = + 1-18N5 
Noy = 5 59Noos Boo = 0.0505No5 
Txx1 = 1-24o. - 1890Ng0Ky, TxxQ 7 0-612Q95 - 66.8 aoxy 


Tyyl = Tyy2 = L.9Roo 


i 


= 94887No5 + GASNgoy* Tyo = - 38.7e°Nap + 33. 4Nogy™ 


Ad 


3. Function Space Approximation. Since so few polynomiais 
give exact solutions, we must have recourses to approximation. 
Gne method of obtaining approximate eolutions is by use of 
function space, * A vector, P, in the function space is the 
ordered set of atress functions, oes 7 T.). The scalar 


yy’ ay 
product of two vectors is given by 


> > 
Pep" 


i] 


Jey, Ty’ + Syy Tyy' + 205yTzy' av 


f 


1 (6, ‘Ts, * ®sy Tyy + 20,4! Thy dav (1) 


fi 


the integration being aver the volume of the body. 
The eatreas state cf un exact solution Ls called a natural 


D> 
state and is Gsnoted by &. 


If the boundary conditions are conditionsa on the atreases, 
or if such displacements 48 aré specified ars zero, the problem 
is said to be one with stress bdDoundary conditions. In this cage 

w@ have also the fsllowing notation: 

3, the completely asscclated atate, satiar&es the equations 
of equilibrium, the surface boundary conditions on the 
atres3es, and the internal boundary stress conditiona. 

5p’ {p=1,..0e8}, the homogene sus associated statea, satisfy 

the equations sf equillbrium, ths internal osundary stress 

eonditions, and have zero stresses at points on the boundary 

where the stresses are apecified. 
ap the orthonormal nomegenssua associated atatesa, are sbtained 


by orthsnsrmaiizeing the eset of 8)’. 


fo Trew AEE 


PEROT TON SPE TI eB we ws 


PARLE TRONS apparel BM! 


6 tet ED 


RE DEI cy Oy STS Re PRN ee) 


SALA OLN TTL BNE fs 


2 Abaya 
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34" (q=l,...,9), the complementer, states, satisfy the 
equation of compatibility, the internel boundary dis- 
placement conditions, and give zeros disvlecements at 
p3ints on the bounccry where zcro Cisplicoments cre 
speciried, 

ty" the orthonormal complementary states, aie abtaines by 

orthonormalizing tue set or ae 

he homogeneous associatce stceres form & Linear mentirsld ™ 
called the h3mx¢enedus essoclited menifold, If eny homo¢gencous 
asa vciated menifoid te acded to the completely assxclated rtate, 
the result ia a completely assrciated state. Th. suvapvace of 
completely Gasxclated states, Mw, is celled the com-letely 
anesocléted subspnee. The complemen&cry states form a lineer 
manifold, Ni, calle? the complementary manifold. M* and N are 
mutually ortho <onal. avery sphere which hos rar ite dlameter 

a line joining « point of Te to oe point or N passes thyowgh 3. 

me uncertsinity in the oaraition ws 5 will be 


y 
least if we pick & point or M* and of such thet the cistance 


between them is €& mini:num. These orints cre «elven by 


> — > o > > ‘ > ; 
Vimn* = 3* - oni! 8" ip )Ip 
> > > > 
Yn" = gay (28 tq")Tq". (2) 


The m&ximunm prssitle error will ve least if we cnr.5e 


G5 an approxinetion for s, the center wf the sphere, friven Ly 


a3 


Aha @ completely associated estate we can take 


Bez (Ta, Tyy*s Txy*) = (2°R0/dy®, a? pe/an®, ~ 2ew/dx2y) (4) 


where 


re 


R(F,*(21) + Fo*(2)) 


n+e 
F\* x malta + iM ote 


= [Eee 
the P,*, Qp*, My*, and Ny,* satisfying the equalt tone (2-2) te (2-5). 
Similerly for our first complementary qia@&d we can take 


> 1 ss ind kb 
By" STs s Tigy fies’) = (0?P,"/dy®, 349," ae, - o°F,"/dx2y) (BEE. 


where F," is defined in 4 manner analsgous tpi Fe e@xeept thet 
equations (2-6) to (2-11) aust be satisfied. 
Then 


8 “8, "7g Wem Ye | (63 
L am rf ae. « 


Using equatione (1-12) and (1-15) we obtain, 


> ‘ t , 
B°-8) = s(52¢aiR] + STSaRTRD + SSS ROR + S3<2RQRQ + SLRTR} 


+ Sok*ka + SPQR, + SoRSRD + 2STR [Ee] + Oy SoRTED 


+ 5 Safe] + OBege2) dev n) 
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Using (1-15), (1-17), and (1-7), equation (7) e&n be written 


t, 


(10) 


P, & " " : " i r) ” 6% 
S°+S, = (2v[Sf( TR +E]E]) + SB(RBRo+ESES) + (Srpr)( HLRZ+RSR)) 
+ 6 85( BT EstESEy) dav (8) 
In a similar manner we sptain 
Beg" = sov[se (attest?) + s2cn8%+232) + a(Seaei RES + 28;6ar "ES ]av (9) 
ang 
> y " 
p°8y = fav [sP( Ry Peel) + 85CRo*e5") + 2Erd RRS + 25,854 Folav 
If we take “? + aM aga our approximation for S, the 
Py mn 2 Pe] s 
maximum error 1a +f" 3". (3 ° 31 )°/87,°3,]*. 


Carrying out the integration »ver the rectangle deseribed 


on page 4, we get 


(-2) 1dPT84e( x, 8; p,q) wkrt 
am + * # " ac Vatia * 
Seah caw ees * S83 (Hep, Hee i 

on s@ . ome: 1s we oo co % aft X 
"Mesp,irteaya? © 4a Pra esp 1 teg, 1° P teg, eed! 

r o ee 7 8 
t Ay (P24) (Tee, Weeg 1% teq, 17 cep, 1) 4 

+ + t " ia 

# Vo(nP*U*t_1) [EP59( GE. oGteg, otTaep, 27 t+q, 2) 


m §B59 (Heep, Hea, 2 +Uisy lteq. 2) + Llp, Q)(Gkey, Herg,2 


rT) : " ) 
*Ht ig oGsep,2) t AolO.9)(Teep erg, otUteg, atcep, dpe (12) 


We heave used — Fo te indigate that ln the summation kK and t 


are to have ths game parity, that le, they ere vath to be sdd 3r 


7 


both even, Similarly for p and q. Tne other new symbols used 


nave the foilowing seaninges: 


15 


- 
» 
~~ 
es 
We’ 
~~ 
i 


57,53) 


For p and q odd 


( Ag(p.a) = SP SAT Ay Jay + TAA) 
| ( PyCp.a) = SPT shit 4s oy + T/A) 
For p and q even ( 
( Ag(p,a) = 8P,83, 


ij 


8(«,t;p,q) 8(p,q;K,t) = O(t,m5q,p) = O(q,p;t, K) 


(krp+2)(trq +2) Ktetl trarh Guay 
Ketel Brqtl « t 
= gt = 4 ‘ ae | 
Gi =e Pi i cl Q13 Tig x $146 My Uys - Sp je%hy 5 


i] 


With the appropriate superacript, * or . In addition, P* and Q* 
are to satisfy equations (2-2) and (2-3), Mm” and N* equations 
(2-4) and (2-5), P" and Q" squations (2-6) to (2-8), and HM" and N" 
equations (2-9) to (2-11). The expressions for 3” B* and 313s 


; ae Beene . u * ry 
ean be sbtained from that 3f S°°S, by repiacing by and by 


ae 


respectively. 
j os “ : ¢ ot 2 : 
Any set sf Fy euch that Trey’ = 877/897 Boy! = r,t /ax ’ 
hs « o°F 0! /2ndy ejue-1 zero om the boundary, end sueh 
' gre cvu.'inervue at y=ere wikl; eerve to 


> 
determine the 3p’. Such @ set 1® given by 


Fp' = (x@ - 0%)?(y? ~DMp*ete.y) £243 
where f(x,y) 19 eoritinusta and Kad aontinusis Terivatives in 
the closed region secupled by the body. 


> 
& set of G 


e 
q 
ments which tace mn zéro values &t points of the boundery where 


Gin be Gonstructes from any sat ow displace- 


the diaplscementsa are specified ts be zero, and which are con- 


tinu sus. 
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It is not necesssry that the F* used to determine the ae 


> ed 
completely agasciated state, S*, satisfy the equation (1-8). ‘f 
any F* such that 1,,% = a“Fe/dy*, Tyy* = 2? pe/ax?, and : 
Tay* = - 2° F#/d x dy satiefy the surface brundary conditions and 4 


a 
guch that "yy 


to obtain an approximate solution. 


* and Thy* are continusus at y = + c san tbe used 


Ag an @xagole, consider the rectangle deseribed on page 4, 


Sudject to the following boundary er 2ndlitiongss: | 


~ 
ret 
rad 
et 
ot 
i 
+ 
is 
% 
tum, 
ey 


(a) te Tea7ey = Oy 


Conditione I are setisfied by taking m = 3 and putting; 


® _ 3% * ‘*. « nt. 5 ot. ~ at. | nt * 
G32 = Hea = Goa = Hao = Gia = Bjo = Taa = Uao = Too = Uge = 0. 
Hoo = -(8/4 + Go), Tho = -Us2 = 0/80n(52, - 585), U2 = wf/8n - The 


From equations (2-2) to (2-5) we get: 


: * ® a8 
5, * 85, > OF) = BE, = 91, = Bh = Th = Phy Me * 8G 0, 


H3) = -(e/4 + G8), TS, = -US, = 0/80n°(S?, - 58), ut, = w/on ~ TH... 


Conditions [1 are eatiefied by 


E . <a ce f ¢ 2 2 
a Paes VG) © “Evay/ S43 .*'2 ay7 


. — — Ei ok wad tka x é “a 
- ? DNS 6 EEE ee POS Tae : os Sg 
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2 
ba 
ho 

i 


4On?(n> - 1) Bn 16(n> - 1) 


K(525 + 535)(n> - 1) «53 
2 l2 22 rad 22 sh = 53.) 


yas aid 


él 


K(S2, + 88,) 68 
7 ? td 6£, . 
hi | rons a =|, (691 - 5) 


| We choose a first aomelemantary atate whieh satisfies the 


surface boundary stresa eonditisna I. This gives; 


| Gbo = R8o 


# 


G22 = HZ2 = Gia = fe = Tho = UB = TE = Ugo = 0. - “4 
Gho = Hfg = -#/8, Tks = -U$p = &/BOn NES - 680), Ula = &/Bn - Tho 
From equations (2-6) to (2-11) we get; 


" — " — tril at vu 
He, = San > Boy = Oly * By > By = Uo, = See Bey =o, 


ros] 

A 

- 
is 


n - © SolSr2 * €a2) - (Sy + Spaéay 
& $4,413 * £o,) - (Syy + So.) ay 


v4 


HY. = & 533 '\2 _* €99) - - (Sy2 + *22)€)) 


8 $41(€y3 + So.) - (St $31) 4.1 


TH) =z ~U%) = ka>/ 80a, n° (6 oP - 63.) 
nf = —* — ("ian 7afen , fat - 9er = Sent] 
S11 - Fay By, eon 
: ee y = Ley. oy 2 3 
ee x ("3 Si2 ~ “25 , fet der aa) ; 
eee ee eee i, Kool $43-5a1)- S2pl€1 yo, )emelerSer-Srta) | . 
te 8Tp-8B 5 ead oat 


&($,.-5,, }-r(é + -€,,5 
, 0844 - Spy) rly y- Say) 184 1509-94 a) 


2hn- 


ee ee ee Le Se Se 
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d= [4,(214-4,)-2)(2%-4p) /ty, 8 = Bolt? A) -) (+H IA, 


re fe A (11-8) -af6(Va-62)Mpag - [AQ¥y- ¥2]/20, 
$45 = Bos y-MoSa gy, Px = Py + 1/4. 
Subetituting these results in equation (11) giver: 

o ay,  {52x08108; + ByuBaHsy* A (OG HB, +HS 68, 143 (Stew ttt 
+ (68, ow* UT UIA tw Naf ur sot et 20468 -w (7h Ty 
m8, Tii) #3(53 =" 9) (UZ, US, US, UE, )+ [5( P,S 7) -w") -5w*( P, -8F3)] 
(aS .0h ULM )* BUS 2i-* *yesa?(P, 88 Vert yus, v8.7} 
r20{ {5/7) (68, +0) +062, (82, +a?) ]88,5,+ ri 7115 $,0" 6) 
a *551 (5510 abo ten isbpdianthomtehshlAy 
~w2 (py om®) (62) 685) ) J 18,05 US iT} * 7 {(n-1) [SF 2@82%$0 
* BBB GalGo* Pol G5 HG2 HG 2052) 195 { [(n-1) Sf o¢(n-1) 97 MeTle 
+ [2 -1)5 3p ¢¢m-1 w2JUt oUt 3+ flin7-1) Pye ( naa )w? J afUTtUtetya)} 
~2 {3 [( 05-1) 6}o-( n-1)w4] (Te oho +Tg ott 5) +5 ((n5-1) 3-(n-1) 94] 
(UP QUgo+UZ aU io) + {5 [(n>-2) Pa8So-(n-1)w 7] -w?(n7-2) (Po+63—)) 
(TPQUSotUS aT} a) + (5 [( 07-2) PoSfo-(n-1)w 4] -5w?( 09-1) (Po-825)} 
(T3aUT +UF a7Sa)} +20/ {(5/7) [in7-1)88+( n-1)wOJew? [in5-1)8 45 
+(n 7 =1)8¢n7}} 18 185° (5/7) ((n? -2)85,+(n-1)w ew? {(n>-1)53, 
+( 9-1) 53007 4}US Ugo (15/7) [(n7-1) 529835 Ayr(ne1)w°} 
#3w° [(n?-1) 8f983o% (n?-1) Pow? ] J-w* [(n9-1) o%(n7 -1)"*}is25+830)) 


(T32U52*U5 2752 }} . 


‘gp 
> 


—= or Ma a 


ig 
Similar re@ults ean be obtained from this equation by substituting 
* for ", or by substituting " for *. 


For spruce with ply of equal thicstness (n=3), the constants 

on page § give: 
won NE 3 — cane A tnt OO EIEN CEFN LIE DAN * a ~ 

4 o/k = (2.26 ~ 9.0435" 4(y/b) - 3.64(y"/b9) + 0.780(x2y/b>) 


Tyxitt = (0.295 - 1.13w*) (y/b) ~ 4, &3(y>/b>) + 0.750( x2y/b3} 
Txx2/t = 3.64 - 10.2(¥/b) - 3.6(99/b3) + 0.750(x2y/d3) 


Txx1/* 


u 


108 ~ 683(y/bj - 125(y3/b2) + 21.2(2%y/b>) 


’ 
Tyyal® = Tyyi/® = Tyyo/k = -0.500 - 0.750(y/b) + 0.250( ¥3/b3) 


Tyyy/% = 2.10 = 17.1 5/d) + 7.05(y9/d3) 
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= 13,,/k = Tryo/® = 0.750(x/b) - 0.750( xy?/b3) 


17.1(x/b) = 21.2(xy*/b>) 


ays 
tq 
™~ 
% 
i 


= (-0.835 + 3.37" - 0,00166w4)/(3070 + 7.9? + 0,00a99"4) 


fs 
{ 


In particular, fer w=30 we get: 
Tyxo/ = 0-406 - 19.0(y/b) - 2.22(y3/b9) + 45.8(x2/a) (¥/d) 
Tyyi/% = 12.1 - 582(y/b) - 16.1(y3/b>) + 274(x2/a%)(y/d) 


Tyyo/% = -0.501 ~ 0.469(y/b) * 0.153(y5/b3) 
Tyyi/® = -0.0153 = 2.18(y/b) + 0.912( 73/3) 


Txy2/* = 4,.59(x/a) - 4.95(x/ally2/b2) 


Tgyl/t = 22.8(x/a) ~ 27.4(x/a)(y?/b?) | 
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4, Approximatiion by use of finite vations. Another 
method of approximmtiion ia to replace the differential equation 
by a get of finite equatisne. To do this we first approximate 
F(x,y) by a polynomial which takea on the values Fay = F(24594) 
at the points (5554) &a shown in (1) below. This can be done 


by uge of Lagrange’s formula of interpolation. 


(1) 


Fo ,i-2 


The etresces are then obtained from (1-2) and from them 
the etraing are »btained by Hopke' law (1-4). These are 


@ubstituted in the compati Pll by equation 


d®exx/ay? & dey, /ax® = 2( 8° ax y/dxdy). (2) 
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This gives 


(-244% 4 + 960, + GOOQ, + 960, - L44a, - 18004, + 8047p 

* Ab2v_) + 452% )Fo g + (9604 - 576%q - 1927, - 1924) F_1 

> (1040_, - 320%, - 480, * 64eL, : 56a ry 2848 5 wn 38485 ~ 3207, 
- 12H IFO * (9604, - 576v_ - 192%, - 1924 )Ps 4 + (86%, 2 of 
+ G44) - 4809, - B24, + LOK, + Z84R) + 3348) - 192H_, eo 
- 320 4 )Fo 4 - (-354 40 + 16024 + POW IF oy + (-38445 

+ 160¥, + 96. )F) 1) + (-384g, * 967, * 160% )F) 3 + (=3844, 


* 96v_) + 160%4)F_1 4 + (-608, + Lddyy - 24v,) - 2h) F_y G 


?- 


(-35a_, * 176a_, + 30a, - 16a, - L19,)Fy _p * (-60,4, + lady) 


- 24y, - 240, Fa 9 * (110% 5 = 16a_, * 30% + 1760, - 350Q)Fo 9 = 0. 


where S% indicates the value of Afor y -.-if y falls on the 
boundary of tro régicns the average of the values 2f @ in tbe 


AX 


two regions--and similarly for #,, 4, and. For a particular 


problem the equations (3) n conven 


@ diagram sigilar to (1) but with the crefficients of Fy y written 
in place of Byoye 
We calculate the value of F and ite first derivatives at 

the boundery by integrating the values of the stresses civen 

on the toundary. The region 1s divided int>s squares of length 
d ty the lines x=x, and Y= 5- The intersections of these 

lines are called nodes. Nodes falling on the boundary, and, 

in cases where the lines Intersect tne boundary in pointe sther 


then nodes, tne nodes jJunt insides the boundary, are called 


boundary nodes. The nodes adjacent to the boundary no@es and 
exterior to the boundary gre Celied adjacent nodes. The values 
of F at adjacent nodes and at boundary nodes which do not fali 
exactly on the boundary are approximated from the values of F 
and ite firat and eecond derivatives at the boundary. Nodes 
interior t32 the boundary which are not boundsry nodes are 
called interior n3deg. An equation can then te sbtsined for 
each interior node by congideriug it as tne point (0, ©) in 
equation (3). Thus the dirferentisl equation has peen replaced 


by the set of finite equations 
«1 oa ; 
ya yF lea yy) oO. (4) 


The solution of these equations can be approximated by 
the eo-called relesstion mothod,© Let F° be the estimated 
value of F at each node (taxing Fo = F for a boundary node 
or an ag@jacent node), Substituting these estimated for F in 


the left side of (4) gives 
J ar yF (x49 9) = RM. (5) 


If KEL = © for all « and l, then F? = F and the estimated 
golution ie the correct eslution. &£f this is not the case, 


ist 


F' (dyes V7) = F°(x~.92) = pad (6 


ad 


The «i'th esuation then becomes 


fi > , «i 9 Ki =, i ] co? a 
BrP (ego) + eee (xysyy) SRE - REY SLR (7) 
jAL i 
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The rest of the equations which invslve F(x;,y¥1) become 
acRe' (xg y1) * Beet FEO (agg y) aR ™ te eXfr' 2 /akt, (8) ad 
JAL 


Thue by shanging the estimated value of F(X, Yz) from F?(x4.¥}) 


pans 


to F'(xe59,), the "residual", Ra has been reduced to ne, 
Starting with the largest IRI, each 18 relaxed tn this 
panner until the error is negligible. The resulting f'a | 


Gi... the approximate adlution. 
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